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Abstract. We study the 'asymmetric' Hubbard model, where hoppings of elec- 
trons depend on their spin. For strong interactions and sufficiently asymmetric 
hoppings, it is proved that the ground state displays phase separation away from 
half-filling. This extends a recent result obtained with Freericks and Lieb for the 
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' Falicov-Kimball model. It is based on estimates for the sum of lowest eigenvalues 
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of the discrete Laplacian in arbitrary domains. 
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[ 1. Introduction 

Electronic properties of condensed matter are difficult to apprehend because of the 
many-body interactions between quantum particles. It is necessary to consider simplified 
£T) ' models that capture the physics of various systems. Of great relevance is the Hubbard 

model |13j where spin-| electrons move on a lattice and interact via a local Coulomb 
repulsion. Although a considerable simplification to the original problem, the Hubbard 
Qh' model is still difficult to study and Hubbard himself considered an approximation where 

^ fH particles of one spin are infinitely massive and behave classically. 

The latter model was reinvented later by Falicov and Kimball in a different context, 
namely in the study of the metal-semiconductor transition in rare-earth materials ,5] . Two 
species of electrons corresponding to different electronic bands are moving on a lattice, 
and relevant interactions are between particles of different species. Electrons carry spins 
^ . but these turn out to be mathematically irrelevant and they can be left aside. There 

exist many results for the Falicov-Kimball model. Let us mention proofs of long-range 
order ^1 ^ H21 EH HB US] and of phase separation ^B] ; all these results are valid 
at half- filling, that is, the total density is equal to 1. Interfaces were studied in [3J. The 
ground state is segregated away from half- filling [01 Hj ; see also [0] . There exist reviews 
by Gruber and Macris J§drzejewski and Lemariski and Freericks and Zlatic [H]. 
Less is known rigorously about the Hubbard model, see the survey by Lieb |22| . 

We consider here a Hamiltonian that interpolates between Hubbard and Falicov-Kimball 
and that describes two species of spinless fermions moving on Z rf ; particles have different 
effective masses, and there is a local interaction involving particles of different species. 
The Hamiltonian in second quantization is 

H A = - c li c yi-t c] x2 c y2 + U^2n x in x2 . (1.1) 

x,ydA x,ydh x£A 

\x-y\=l \x-y\=l 

Here A is a finite cube in Z d and c\- and c x j are creation and annihilation operators of 
a fermion of species j at site x. The first two terms represent the kinetic energy of light 
and heavy electrons respectively (we suppose that ^ t ^ 1). n x j = c\-c x j is a particle 
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number operator. The positive parameter U measures the strength of the on-site repulsion 
between particles of species 1 and 2. 

Setting t = 1 yields the Hubbard model, and t = yields the Falicov-Kimball model. It 
is interesting to note that the behavior of both models is similar when both particles have 
density ^: for d ^ 2, the ground state of the Hubbard model is a spin singlet and 
the one of the Falicov-Kimball displays long-range order of the chessboard type ^7j. This 
holds for all strictly positive values of the coupling constant U. It is natural to conjecture 
that long-range order occurs for all t. 

Convergent perturbative expansions for large U are a major source of results for the 
Falicov-Kimball model, at least at half-filling. See [3] and ^H] for general methods, and 
[2] for a discussion specifically to the Falicov-Kimball model. These methods are robust 
and extend to any perturbation of the model. This holds in particular in the case of the 
asymmetric Hubbard model with small t. 

Our goal is to identify a phase with with segregation and to contrast it with chessboard 
order and with high-temperature disorder. This suggests to look at the following operator, 

O/Jp) = yn ~ n V+xd 2 - ( L2 ) 

1 1 yeA 

The corresponding correlation function is given by the expectation of o~a(x) in the equi- 
librium state. We consider here the canonical ensemble where densities of light and heavy 
particles are fixed to p\ and p2 respectively. High temperature states are translation in- 
variant and exponentially clustering, and the correlation function converges to 2/92(1 — P2) 
as x — > 00. Notice that < 2/92(1 — P2) ^ \- We identify here a domain of parameters 
where the expectation of cr A (x) is zero in the ground state (segregation). At half-filling 
perturbation methods [HI Hj show that it is close to 1 when \x\ is odd (chessboard 
order). 

Theorem 1. Suppose that p\ + P2 ^ 1- There exist Uq < 00 and to > (that depend on 
pi and p2 only) such that for U > Uq and t < to we have 

(T A ,a A (x)T A ) =0(|A|-*). 

Here T\ is any ground state in the subspace where light and heavy particles have densities 
pi and p2, respectively. 

This theorem extends the result of [SI E| for the Falicov-Kimball model. Its proof pro- 
ceeds by obtaining estimates for the ground state energy. The ground state is a linear 
combination of states with a fixed configuration of heavy particles. The weight of config- 
urations with large 'boundary' (pairs of nearest-neighbor sites where one is occupied and 
one is empty) is small. Indeed, most of light particles are delocalized in the remaining 
sites and their kinetic energy would otherwise be great, as it is roughly proportional to 
the boundary (see and Section [2j) . The pressure exerted by the light particles packs 
the heavy particles together. The kinetic energy of heavy particles is therefore irrelevant, 
and simple estimates suffice in bounding their contribution. These ideas are detailed in 
Section 0J 

Section 12 reviews the results for the sum of lowest eigenvalues of the discrete Laplacian 
obtained in with some improvements in the regime of low densities. We discuss the 
segregated states of the asymmetric Hubbard model for all ^ t ^ 1 in Section |3 For 
given densities of light and heavy particles, there is one free parameter to characterize 
segregation: the proportion of volume occupied by each type of particles. The restricted 
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phase diagram of segregate states displays a transition between a phase where the local 
density of heavy particles is maximum (that is, 1), and a phase where they have a local 
density that is strictly less than 1. Section 0] is devoted to the proof of Theorem ^ 

Acknowledgements: I am grateful to both referees for useful comments, and to Loth 
Hermi for drawing my attention to the paper |24| . This paper is dedicated to Elliott Lieb 
on the occasion of his seventieth birthday. With Tom Kennedy, Elliott Lieb reinvented the 
Falicov-Kimball model and obtained the first rigorous results ^7]. I enormously benefitted 
from a collaboration with him and Jim Freericks on segregation in this model; the present 
paper is directly inspired by [5]. 



2. Sum of lowest eigenvalues of the discrete Laplacian 

The sum of the iV lowest eigenvalues of the discrete Laplacian in a finite domain AcZ rf 
gives the ground state energy of N spinless, non-interacting electrons hopping in A. This 
quantity is relevant to some problems of condensed matter physics. Thermodynamics 
suggests that it is equal to a bulk term that is proportional to the volume |A| of the 
domain, plus a positive boundary correction that is proportional to the boundary of A. Li 
and Yau proved in 1983 that the sum of lowest eigenvalues of the continuum Laplacian is 
indeed bounded below by the bulk term [20] • See [231; Theorem 12.3, for a clear exposition. 
The proof readily adapts to the case of the lattice. 

The problem on the lattice turns out to be simpler and allows for bounds on the bound- 
ary correction, for given 'electronic density' p = N/\A\. Precisely, the boundary correction 
can be bounded above and below by positive numbers times the 'surface' of the boundary. 
This was done in .6]; this section contains some improvements in the limit of low densities. 

Corresponding statements in the continuum case have not been obtained yet. The best 
statements seem to be the upper bound of Lieb and Loss, Theorem 12.11 in [22], and the 
lower bound of Melas [23], who obtained a positive correction of the order of the size of the 
domain to the power d — 2. However, these bounds are not proportional to the boundary 
when the density is fixed. 

For a finite domain AcZ d the discrete Laplacian h\ is defined by 

h A ip( x ) = - f(y) + 2 Mx), (2.i) 

\y-x\=l 

for all x G A. Here (p G li (A) is a normalized, complex function on A. If <p is an eigenstate 
with eigenvalue e, so is (— l)^ip with eigenvalue Ad — e (here \x\ denotes the l\ norm of 
x G Z d ). One also checks that h\ #J 0, and therefore its spectrum is contained in (0,4d) 
and is symmetric around 2d. The bulk term involves the ground state energy per site e(p) 
of free fermions and it is given by 

<P) = 77^77 f £ kdk. (2.2) 



{^Y Je k <e F (p) 

Here k G [-7r,vr] d and e k = 2d - 2 £? =1 cos k, L . Notice that |A;| 2 - ^|A;| 4 < e k < \k\ 2 . The 
'Fermi level' £f(p) is defined by the equation 

P=7^f dk - (2-3) 



(2vr) c 



e k <e F (p) 
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Let -B(A) denote the number of bonds connecting A with its complement, 

B(A) = \{(x,y) : x G A,y G A c , and \x - y \ = l}|. (2.4) 

If £A,iV is the sum of the N lowest eigenvalues of h A , and p = N/\A\ is the density, we are 
looking for bounds of the form 

e(p)\A\ + a(p)B(A) < S AiN < e(p)|A| + b(p)B(A) (2.5) 

with positive a(p), 6(p), that are independent of the domain. It was proved in |H| that 

Kp)=P~tAp) ( 2 - 6 ) 

gives the optimal upper bound, that is saturated by domains consisting of isolated sites. 
(The size of the boundary was defined differently in jS] but minor changes in the proof 
yield the upper bound stated here.) 

We define a(p) to be the minimal 'surface energy' among all possible domains. Namely, 
for p € [0, 1] n Q, 

/ n ■ r S A ,N-e(p)\A\ 
a( P )= mf m • (2.7) 

The infimum is taken over all finite domains A such that p\A\ = N is an integer. The 
symmetry of the spectrum of h A around 2d implies that a(l — p) = a(p). We give below 
lower and upper bounds, stating in particular that a(p) > for < p < 1. Many questions 
remain open, such as the existence of a minimizer in (|2.7jl ; continuity of a(p); monotonicity 
and convexity of a(p) for ^ p ^ ^. It is even not clear whether the infimum (|2.7|) can 
be taken on connected sets. In order to state the bounds for a(p), let us introduce 

Z(j>) = pe F (p)-e(j>). (2.8) 

Theorem 2. For all < p ^ ^, we have 

< a{p) ^ l d i(p). 

For small densities, we have 

a(p)>j£ F £(p)(l-0(p 2 / d )). 

We prove here that a(p) is bounded below by j^jC(p) a t low densities and that it is 

smaller than -^C(p)', notice that £(p) ~ p 1+ ~a as p — > 0. Efforts are made here to get the 
best possible factor. On the other hand, pushing the range of densities instead, we could 
get a positive lower bound for < p < (4-7r)~ d/ ' 2 r(| + 1) -1 . Remaining densities are much 
more difficult to treat and we refer to [£] (and to |Hj for subsequent improvements and 
simplifications) . 

Proof of the lower bound for a(p) for low densities. We follow <o , with some improvements. 
Let ipj be the eigenvector of h A corresponding to the j-th eigenvalue ej, and tpj be its 
Fourier transform 

<Pj(k) = \>,0r)e ifca \ ke[-7r,7r] d . (2.9) 



icgA 



Then 



Sa ' N = J2^I .P( k ) £ kdk, (2.10) 
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where 

N |A| 

p(k) = j2\^)\ 2 = \M- E i&(*)i 2 - (2- 11 ) 

3=1 j=Af+l 

We also observe that , 2 J p(k)dk = N. One obtains a lower bound for £a,at by taking 

the infimum of the right side of Q2.1UJI over all positive functions p smaller than |A| 
and with the proper normalization. This gives the bulk term |2l)l 1231 E|. In order to 
extract the effect of the boundary, one strengthens the upper bound for p(k), aiming at 
|A| — const • .B(A). We start as in |H] and write down a Schrodinger equation that is valid 
for all x € TL d : 

— tfj (x + e) + 2d<fj (x) + Xac (a) ^ <Pj( x + e ) = e j<Pj( x )- ( 2 -l 2 ) 

e e:a;+egA 

It is understood that <fj(x) = if x ^ A; the sums are over unit vectors e. The term with 
the characteristic function Xac involves only sites that are close to the boundary. The 
Fourier transform of this equation can be written as 

e k ifj(k) + (b k ,<pj) = ej0j(k), (2.13) 

where bk is the following 'boundary vector' 

b k (x) = X 9A (x)e- ihx £ e" ifce . (2.14) 

e:x+efA 

We introduced the set dA of sites inside A touching its complement 

<9A = {x € A : dist (x, A c ) = 1}. (2.15) 

We observe that B(A) ^ ||6fc|| 2 ^ 2dB(A), the lower bound holding at least when |fc|oo ^ j. 
The last term of (|2.11|) can then be written using (|2.13|) as 

j=Ar+i j=Ar+i £fc J (EiJjv+1 l e i - £ fcl 2 l(Wj)l 2 ) 

The lower bound follows from Holder's inequality. One easily checks that 

|A| 

E \ e i _£fc l K^fe'^i)! 2 ^ (h,hAb k ) - (e k + ejv)||& fc || 2 , 

(2.17) 

|A| V ! 

E \ e j ~ £ k\ 2 \(b k ,ipj)\ 2 < (b k ,h 2 A b k ) + e\ \\b k \\ 2 . 
j=N+l 

From now on we suppose e& and ejy to be small so that they add to less than 1. Notice 
that (b k ,h\b k ) = \\h\b k \\ 2 . Because each site of dA has a neighbor outside A and b k is 
zero there, we have 

(b k , h k b k ) = £ \b k {x) - b k {y)\ 2 > £ |6 fc (x)| 2 = ||6 fc || 2 . (2.18) 

{a;,i/}:|x— j/|=l x£dA 

Then (|2~TT)|) . (|2~T7j) . and (t2~THl) . imply that 

El» n-M2 ^ (b k ,h A b k ) 4 (l - e k - cat) 4 
* (l|/uM 2 + . 2 |N| 2 )3 • (2 ' 19) 
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We estimate the denominator. 

IIM>*II 2 = E| E [h(x)-b k (y))\ 

x&A y£Z d 
\y-x\=l 

^ 2d E E \ b k(x)-h(y) 

xeA yG z d 
\y-x\=l 



(2.20) 



2d E \h(x)-h(y) -2d \ b k(y)f 



x,yeZ d 
\x-y\=l 

^4d(b k ,h A b k )-2d\\b k \\ 2 . 
Inserting this bound in l|2.19|) . we obtain 

|A| 

E \<Pj( k )\ 2 > i 1 - ?k - e N ) 4 \\b k \ 



xeA c ,y£A 
\x~y\=l 



(b k ,hAh)/\\b k \ 



j=N+l 



[Ad - (2d 



-2^ WP 13' 

''h){b k ,h A b k )l 



(2.21) 



Simple analysis shows that the minimum of the fraction under the condition (b k , h\b k ) ^ ||frfc|| 2 
\h(x)\ 2 



is equal to 73^3- (1 — f|)- Furthermore, b k (x) is close to bo(x) for small k, 



(cos ke — i sin ke) ^ | cos ke 

e:x+e^A e:x+e<£A 
2 



>|(l-|e*) E x | X 1 

e:x+e4A 



(2.22) 



Clearly, (1 - e k - e N ) 4 (l - ||)(1 - e k ) > 1 - 6e fe - 4 eA r. Using (t2~TH and lt2~32l and since 



|6o|| 2 > 5(A), we get 



|A| 

E \<Pi&)\ 2 > j^(l-^ k -4e N )B(A). 

j=N+l 



(2.23) 



We can insert this estimate into 1)2.1 ljl so as to get 

p(k) < |A| - (^(1 - 6e k - 4e N )B(A). (2.24) 

Suppose we have a bound p(k) ^ (1 — a)|A| for some a that is independent of k. Lieb 
and Loss 'bathtub principle' (Theorem 1.14 in [23]) yields 



S W > (1 - «)l A l(^p jf 

= (l-«)e( T ^)|A 



£fc d/c 



(2.25) 



And because (1 — a)e(j3^) is convex as a function of a, and that its derivative is equal 
to £(^7-), we obtain 

S KN ^e(p)\K\ + a^{p)\K\. (2.26) 

Let < rj < 1. We define 



(33J 



(2.27) 
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The condition (|2.24|) implies that p(k) ^ (1 - a)|A| for all k such that e k < e F (j^), 
provided the following condition holds true, 

??^6e F ( T ^) + 4e 7 v. (2.28) 
Given rj, we restrict to densities small enough so that 

^{jz E ^)<h- ( 2 - 29 ) 

27d 2 

(Notice that is an upper bound for a.) Then for all domains A and all numbers of 
electrons N such that 4e^ < |r/, the condition (|2.28|) is satisfied and we obtain 

S AtN > e(p)\A\ + (^(1 - mp)B(A). (2.30) 

Consider now the case where (|2.29j) is fulfilled but Aejq > \r\. We define N' such that 
4ejv ^ \r\ and 4ejv'+i > \r\. Then 

N 

S\,n = S\,N' + ej 

j=AT'+l (2-31) 

> e(p')|A| + (3^(1 - riWWA) + Hp - p')\K\. 
The right side is larger than e(yo)|A| + 73^(1 — v)C{p)B(A) provided that 

e(p') + (3^(1 - ^^pOiaT + HP - P') > <P) + w(! " VMP)^- (2-32) 
It is enough to check that the function ~rjp — e(p) — 2fdF^(p) * s increasing. The derivative 
of £(p) is equal to p^ey(p)- It is possible to verify that 

i- p e F (p) < f T(| + lf' d p- 1+1 d (2.33) 

(the bound is optimal in the limit p — > 0). The function above is therefore increasing for 
small densities. The number rj can be chosen arbitrarily small by taking the density small 
enough. Precisely, the condition is that const • p 2 / d ^ r/. This means that given p, we can 
take 7] = 0{p 2 / d ). □ 

Proof of the upper bound for a(p). Let A be a (rather large) domain, and A' be a set of 
isolated sites outside of A. Let r be such that [A'| = r|A|. The spectrum of /iauA' 
is given by the union of the spectrum of h\ and of {2d}, the latter eigenvalue being 
at least |A'| times degenerated. We have 5a,jv ^ <SauA',jv (with equality if M ^ |) and 
B(A') = -B(A) + 2dr|A|. Using the upper bound for Sa.a and the lower bound for 5auA',a 5 
we obtain (with p = | A | (^ +r ) ) 

e((l + r)p) |A| + 6((1 + r)p)B(A) > (1 + r)e(p)|A| + a(p) [5(A) + 2dr|A|] . (2.34) 
Reorganizing, 

a(p) [^p + 2dr] < e((l + r)p) - (1 + r) e(p) + 6((1 + r)p)^-. (2.35) 

This inequality holds for any domain A such that (1 + r)p|A| is an integer. Ratios bound- 
ary/volume can be made arbitrarily small and therefore the corresponding terms can be 
omitted. We obtain 

e((l + r)p) -(l + r )e{p) 
a{p) < — . (2.36) 

Taking the limit r — > yields the result. □ 
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3. A DISCUSSION OF SEGREGATION 

Particles of different species segregate away from half-filling, at least for large U and 
small t. The domain splits into two subdomains, A = Ai U A2, with Ai containing light 
particles only, and A2 containing heavy particles only. This is true up to boundary terms 
that do not contribute to the bulk energy. We neglect boundary terms in this section. 

There is one free parameter that controls segregation, namely the ratio of the volumes 
occupied by each phase. For t = the ground state is realized with IA2I = N%; it was 
argued in [Hj that light particles exert a 'pressure' that packs heavy particles together; this 
pressure overcomes the tendency of heavy particles to delocalize so as to decrease their 
own kinetic energy. If t is large enough however, heavy particles will extend their domain. 
We study this mechanism in this section, assuming that particles always segregate. From 
the point of view of rigorous results, we obtain upper bounds for the ground state energy 
of the system. 



Ai \ 




pi 




l-v 






A 2 




P2 




V 



A 



Figure 1. A segregated state involves a partition of the domain into sub- 
domains Ai and A2 for light and heavy particles respectivey. The boundary 
between subdomains is supposed to be small so that its contribution to the 
bulk energy is negligible. With v such that |A 2 | = ^|A|, densities inside 



each sub domain are and 



We consider a finite domain A C Z rf partitioned in two subdomains Ai and A 2 . We 
fix the number of particles N± and N 2 of light and heavy particles respectively, and we 
denote the corresponding densities by p\ = iVi/]A| and p 2 = AT 2 /|A|. Let v = |A2|/|A|; we 



have P2 ^ v ^ 1 — p\ . Notice that the densities inside each subdomain are = and 
. Neglecting the contribution of boundaries, the energy per site of this segregated 



N 2 _ P2 
[Ai] V 

state is 



e(p 1 ,p 2 ;v) = (l-v)e(&)+tve(%). (3.1) 

For given densities p\ and p 2 we are looking for the minimum of e(p\, p 2 ;i^) with respect 
to v. One easily computes 

^e(p uP2] u)=^)-taf). (3-2) 

It is worth noticing that e(pi, P2',v) is convex in u, as its second derivative is positive (£ 
is increasing). At v = p 2 , we have 



^-e(p 1 ,p 2 ;v) 



V=p2 



Pi 



-pi 



2dt. 



(3.3) 
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This expression can be positive or negative, the critical parameter being t c = ^Ci i- )• 
On the other hand, the derivative of e(pi,p2] v) at v = 1 — p\ is always positive (if t ^ 1). 
Therefore the segregated state that has minimum energy (among segregated states) is 
given as follows: 

• If t ^ t c = ^C(tt^)! the minimizer is v = p2, and the phase of heavy particles 
has density 1. 

• If t c < t < 1, the minimizer v is between p2 and pi P ^_ p2 , and the phase of heavy 
particles has a density strictly larger than p\ + p2 and strictly smaller than 1. 

• If t = 1 the minimizer is v = pi P + p2 and the phase of heavy particles has density 

Pi + P2- 

This is illustrated in Fig. |2 that displays a restricted phase diagram where only segre- 
gated states are considered. This description is relevant only if a segregated state minimizes 
the energy. This is proved in the case of small t. 

A major open question in this model is whether segregation really occurs for t > t c . 




Figure 2. Restricted phase diagram for segregated states. The phase of 
heavy particles has density 1 in the dark gray domain; its density is strictly 
less than 1 in the light gray domain. Segregation is proved in the black 
domain when U = oo (and in a smaller domain when U is large). 



4. The ground state of the asymmetric Hubbard model 

Let .F(A) be the Fock space for spinless fermions in A. For A' C A, let $a' £ -^(A) 
represents the state with |A'| particles occupying all sites of A'. {^a'Ia'cA is a basis for 
JF(A). The state space for the asymmetric Hubbard model is J^(A) (g>.F(A). Any function 
Ta € ^"(A) ® J^(A) can be written as 



Ai,A 2 cA 



(4.1) 
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witn Eai,a 2 l a Ai,A 2 | 2 = !• Let a A 2 = (Eaj I«a 1: a 2 | 2 ) 1/2 and *(A 2 ) G .F(A) be the nor- 
malized function such that 

«A 2 ^(A 2 ) = ^.A^Ai- (4.2) 
AiCA 

Then ^a 2 a\ 2 = 1, and the function Ya can be written as 

T A = «a 2 ^(A 2 )0$a 2 - (4.3) 

A 2 cA 

We derive in Proposition|3]below an inequality for the coefficients a\ 2 that will allow us 
to establish segregation in the ground state of the strongly asymmetric Hubbard model. 

Let J- (A; N) denote the Hilbert subspace of ^"(A) corresponding to particles. That 
is, it is spanned by {^A'} with |A'| = N. All spaces J-(A;Ni) ® J-(A; N2) are invariant 
under the action of H\ since the latter conserves both particle numbers. As before, we 
denote densities by p\ = rA and i»2 = pa- The term 7(C) that appears below was defined 

in [H]; it behaves like ^ for large U. 

Proposition 3. Let T A be a ground state of H A in T(A;Ni) <£> J 7 (A; iVjj). If cl(jz^) > 
j(U) + t, we have 

^ a ^ (A2 ^ a (T ^)-V)-t |A| 
A 2 cA ^1-P2 J lx ' 

Proof. We write the energy of a state using coefficients a\ 2 defined above, and then use 
results obtained for the Falicov-Kimball model. Let Ta be the kinetic energy operator for 
particles in A; it acts on .F(A), and can be written as 

t a = - E c * c v ( 4 - 4 ) 

x,ydA 
\x-y\=l 

where c x and c x are creation and annihilation operators of a fermion at x. Notice that the 
kinetic terms of are given by Ta <8) H + tt ® Ta- Furthermore, for A' C A, let Va A' 
be the operator 

V A , A ' = UJ2 4^- (4.5) 

zeA' 

It represents an external potential that is equal to U on sites of A' and otherwise. The 
energy of a state Ta given by Q4.3JI can be written as 



(? A ,H A ? A )=t £ aA 2 aA' 2 (*(A 2 ),*(A^))($A 2! TA*Ai 



A 2 ,A 2 



^a| 2 (*(A 2 ), (T A + F a ,a 2 )*(A 2 )). (4.6) 



A 



(The sums are over sets satisfying |A 2 | = \A' 2 \ = A r 2 .) Notice that the first term of the 
right side is bounded by 

t £ >a A2 a A , <t(£ '4 2 ) 1/2 ( E '"if' =tJ2 a l B ^ ( 47 ) 



A 2 ,A' 9 A 2,AL A2A9 
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where the symbol means a sum over pairs of sets A 2 ,A 2 C A that differ only by 
one site moved to a neigbor (that is, the symmetric difference A2AA2 must be a pair of 
nearest-neigbors) . 

The strategy of the proof is to consider the expression Q4.6J1 for the energy of the ground 
state Ta- We get an upper bound by using a trial function that is independent of the 
coefficients {a\ 2 }. We then estimate the second term of (j4.6|) from below, using the 
inequality ()2.5|) for the segregation energy. The corresponding expression involves the 
coefficients {a\ 2 } and we obtain the inequality stated in Proposition |3J 

Let A2 C A be such that IA2I = A2, and let us consider ^^2) <8> ^a 2 where ^(A2) is a 
normalized function of ^"(A, Aq) with support on A 2 . We have 

(f(A 2 ) ® <5>a 2 ,H a tf(A 2 ) ® $ Aa ) = (tf (A 2 ), T A c^(A 2 )) 

^|A|| e ( T ^)+6(T^)^(A2). (4.8) 

We used the upper bound in (|2.5jl . We take for A2 a square if possible, or a domain with 

l 1 

very close shape. Its boundary is less than 4dN 2 d ■ The boundary term in (|4.8|) is then 

l—- 1 

smaller than Adb( 1 p } p2 )p 2 d \A\~d. We use now the lower bound in (|2.5f) . As stated in 
this paper it holds only for U = 00. However, it was extended in jS, to finite U; namely, 
it was shown there that 

(*(A 2 ), (T A + Va)*(A 2 )) > e(^)|A§| + [o(^) - i{U)]B{A 2 ), (4.9) 

where o(-) is the minimal surface energy defined in (|2.7[) . Combining this with Q4.6|) and 
1)4.71) . we have for any ground state function Ta, 

(T A ,iJ A T A )^-i^al 2 S(A 2 ) + ^aiJ e (A ) |A 1 | + ( a (_£L_)_ 7(C /)) B( A 2 )]. (4.10) 

A 2 A 2 

We now compare this expression with the upper bound Q4.8JI and we obtain 

J2al 2 B(A 2 )[a(j^- 2 )- 7 (U)-t] < Ad Kj^) pl~ ^ \A\^. (4.11) 

A 2 

□ 

Proof of Theorem^ Using the decomposition (|4.3|) for the ground state T A , we have 

(T A ,a A (x)T A ) = E al 2 X A2 (y)X A c(x + y). (4.12) 

1 y£k A 2 cA 

52x A2 (y)X A c(x + v)^B(A 2 )\x\ o0l (4.13) 

y£A 

(T A ,a A (*)T A ) ^ E 4 2 ^(A 2 ). (4.14) 

' ' A 2 cA 

Hole-particle symmetries in this model are similar to those in the Falicov-Kimball model, 
see ^7j, and allow to restrict to the case p\ + p 2 < 1. We have a( 1 ^ 1 2 ) > 7(C) + i if f/ 
is large and t is small, so that Proposition El is valid. We use it to control the sum above 
and we get Theorem ^ □ 



It is clear that 



and therefore 
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